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The  modeling  of  the  pressure-strain  correlation  of  turbulence  is  examined  from  a  basic 
theoretical  standpoint  with  a  view  toward  developing  improved  second-order  closure  models. 
Invariance  considerations  along  with  elementary  dynamical  systems  theory  are  used  in  the 
analysis  of  the  standard  hierarchy  of  closure  models.  In  these  commonly  used  models,  the 
pressure-strain  correlation  is  assumed  to  be  a  linear  function  of  the  mean  velocity  gradients 
with  coefficients  that  depend  algebraically  on  the  anisotropy  tensor.  It  is  proven  that  for 
plane  homogeneous  turbulent  flows  the  equilibrium  structure  of  this  hierarchy  of  models 
is  encapsulated  by  a  relatively  simple  model  which  is  only  quadratically  nonlinear  in  the 
anisotropy  tensor.  This  new  quadratic  model  -  the  SSG  model  -  is  shown  to  outperform  the 
Launder,  Reece,  and  Rodi  model  (as  well  as  more  recent  models  that  have  a  considerably 
more  complex  nonlinear  structure)  in  a  variety  of  homogeneous  turbulent  flows.  However, 
some  deficiencies  still  remain  for  the  description  of  rotating  turbulent  shear  flows  that  are 
intrinsic  to  this  general  hierarchy  of  models  and,  hence,  cannot  be  overcome  by  the  mere 
introduction  of  more  complex  nonlinearities.  It  is  thus  argued  that  the  recent  trend  of 
adding  substantially  more  complex  nonlinear  terms  containing  the  anisotropy  tensor  may  be 
of  questionable  value  in  the  modeling  of  the  pressure-strain  correlation.  Possible  alternative 
approaches  are  discussed  briefly. 


‘This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NAS1-18605  while  the  first  and  second  authors  were  in  residence  at  the  Institute  for  Computer 
Applications  in  Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 


1.  INTRODUCTION 


The  pressure-strain  correlation  plays  a  pivotal  role  in  determining  the  structure  of  a  wide 
variety  of  turbulent  flows.  Consequently,  the  proper  modeling  of  this  term  is  essential  for 
the  development  of  second-order  closure  models  that  have  reliable  predictive  capabilities. 
Rotta  (1951)  developed  the  first  simple  model  for  the  slow  pressure-strain  correlation  (i.e. , 
the  part  that  is  independent  of  the  mean  velocity  gradients)  which  describes  the  return 
to  isotropy  behavior  of  turbulence  within  the  framework  of  a  full  Reynolds  stress  closure. 
This  model  has  served  as  a  cornerstone  for  the  representation  of  the  slow  pressure-strain 
correlation  in  a  variety  of  the  commonly  used  second-order  closures  such  as  the  Launder, 
Reece,  and  Rodi  (1975)  model.  Subsequent  to  this  work,  Lumley  (1978)  demonstrated  the 
need  for  nonlinear  terms  in  models  for  the  slow  pressure-strain  correlation  and  derived  a 
nonlinear  representation  theorem  for  the  slow  pressure- strain  correlation  based  on  isotropic 
tensor  function  theory.  In  the  high-Reynolds-number  and  small  anisotropy  limit,  the  Lumley 
(1978)  model  reduces  to  the  Rotta  model. 

The  simplest  model  for  the  rapid  pressure- strain  correlation  that  is  used  in  second-order 
closure  modeling  is  based  on  the  assumption  of  isotropy  of  the  coefficients  of  the  mean 
velocity  gradients;  this  gives  rise  to  a  rapid  pressure-strain  model  with  a  single  term  that 
is  proportional  to  the  mean  rate  of  strain  tensor  (see  Rotta  1972  and  Mellor  and  Herring 
1973).  Starting  with  the  work  of  Launder,  Reece,  and  Rodi  (1975),  anisotropic  models  for 
the  rapid  pressure-strain  correlation  have  been  formulated  wherein  the  coefficients  of  the 
mean  velocity  gradients  are  taken  to  be  functions  of  the  anisotropy  tensor.  In  the  Launder, 
Reece,  and  Rodi  model,  the  fourth-rank  tensor  of  coefficients  of  the  mean  velocity  gradient 
tensor  is  linear  in  the  anisotropy  tensor  whereas  most  of  the  newer  models  developed  during 
the  last  decade  are  nonlinear  (see  Shih  and  Lumley  1985,  Haworth  and  Pope  1986,  Speziale 
1987,  Reynolds  1987,  and  Fu,  Launder,  and  Tselepidakis  1987).  The  nonlinear  models  of 
Lumley  and  co-workers  have  been  primarily  developed  by  the  use  of  realizability  constraints 
(see  Lumley  1978).  In  contrast  to  this  approach,  the  nonlinear  model  of  Speziale  (1987)  was 
derived  using  a  geostrophic  flow  constraint  (i.e.,  material  frame-indifference  in  the  limit  of 
two-dimensional  turbulence)  whereas  Reynolds  (1988)  has  attempted  to  develop  models  that 
are  consistent  with  Rapid  Distortion  Theory  (RDT). 

In  this  paper,  the  general  hieiarchy  of  closure  models  for  the  pressure-strain  correlation 
will  be  considered  which  are  linear  in  the  mean  velocity  gradients,  with  coefficients  that  are 
functions  of  the  anisotropy  tensor.  This  hierarchy  of  models,  which  was  motivated  by  analy¬ 
ses  of  homogeneous  turbulence,  encompasses  all  of  the  closure  models  for  the  pressure-strain 
correlation  that  have  been  used  in  connection  with  second-order  closures.  A  general  repre¬ 
sentation  for  this  hierarchy  of  closure  models  will  be  derived  by  means  of  invariant  tensor 
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function  theory.  This  general  representation  for  the  pressure-strain  correlation  will  then  be 
applied  to  plane  homogeneous  turbulent  flows  -  the  class  of  flows  that  have  long  played  a  piv¬ 
otal  role  in  the  screening  and  calibration  of  such  models.  However,  there  will  be  one  notable 
difference  with  previous  work  on  this  subject:  the  simplest  generic  form  of  this  hierarchy  of 
models  will  be  sought  that  has  the  same  equilibrium  structure  in  the  phase  space  of  plane 
homogeneous  turbulent  flows  as  the  general  model.  This  generic  form  -  which  will  be  termed 
the  SSG  model  -  is  only  quadratically  nonlinear  in  the  anisotropy  tensor.  It  has  the  advan¬ 
tage  of  being  topologically  equivalent  to  the  general  model  in  plane  homogeneous  turbulence 
with  the  simplicity  of  structure  that  allows  for  the  determination  of  all  empirical  constants 
based  on  calibrations  with  pertinent  RDT  results  and  two  well  documented  physical  exper¬ 
iments  (i.e.,  homogeneous  turbulent  shear  flow  and  the  decay  of  isotropic  turbulence).  This 
new  SSG  model  will  be  shown  to  perform  better  than  the  commonly  used  Launder,  Reece, 
and  Rodi  model  for  a  variety  of  homogeneous  turbulent  flows  which  include  plane  strain, 
rotating  plane  shear,  and  the  axisymmetric  expansion/contraction.  However,  there  are  still 
some  remaining  deficiencies  in  the  new  model,  particularly  for  rotating  shear  flow.  Based 
on  an  analysis  of  the  bifurcation  diagram  for  rotating  shear  flow,  it  will  be  shown  that  these 
deficiencies  are  intrinsic  to  this  general  hierarchy  of  pressure- strain  models  and  cannot  be 
eliminated  by  the  addition  of  more  complex  nonlinear  terms.  The  implications  that  these 
results  have  on  turbulence  modeling  will  be  discussed  in  detail  along  with  suggested  future 
directions  of  research. 


2.  THE  GENERAL  PRESSURE-STRAIN  MODEL 


We  will  consider  the  turbulent  flow  of  a  viscous,  incompressible  fluid  governed  by  the 
Navier-Stokes  and  continuity  equations 


dv{  dv{  dP  , 

Tt  +  +  ^ 


dvi 

dx{ 


=  0. 


(1) 


(2) 


In  (I)  -  (2),  Vi  is  the  velocity  vector,  P  is  the  modified  pressure,  and  v  is  the  kinematic 
viscosity  of  the  fluid.  The  velocity  and  pressure  are  decomposed  into  ensemble  mean  and 
fluctuating  parts,  respectively,  as  follows: 


V{  =11,-1-11,,  P  =  P  +p. 

Here,  the  mean  and  fluctuating  velocity  are  solutions  of  the  transport  equations 

dvi  _  dvi  dP  _2_  <9t,. 
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where  T,j  =  ~U{Uj  is  the  Reynolds  stress  tensor. 

The  Reynolds  stress  tensor  T;j  is  a  solution  of  the  transport  equation 

dvj 


dr,,  _  drij 
+  vk- 


dt 


'dxk 


-Tik 


dxk 
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Tjkdxk+  dxk  n,J  +  £lJ  +  vV  T,J 


where 
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(9) 

(10) 

(11) 


are,  respectively,  the  third-order  diffusion  correlation,  the  pressure-strain  correlation,  and 
the  dissipation  rate  correlation.  Equation  (8)  is  obtained  by  taking  the  symmetric  part  of 
the  ensemble  mean  of  the  product  of  the  fluctuating  velocity  Uj  with  equation  (6).  For  homo¬ 
geneous  turbulent  flows,  at  high  Reynolds  numbers  where  the  dissipation  is  approximately 
isotropic,  the  Reynolds  stress  transport  equation  (8)  simplifies  to 


dvj  dv{  2 

T"  -  -T*d7k  ~  TikdTt  n”  f  3e6'’ 


(12) 


where 

du,  dui 
dxk  dxk 

is  the  scalar  dissipation  rate.  Equation  (12)  becomes  a  closed  system  for  the  deterrrdnation 
of  T,j  in  terms  of  dVi/dxj  once  closure  models  for  II, j  and  t  are  provided.  Sine-  is  the 
only  unknown  correlation  that  contains  directional  information,  it  is  clear  that  it  will  play 
a  pivotal  role  in  determining  the  structure  of  the  Reynolds  stress  tensor  for  a  given  mean 
velocity  field.  This  dominant  influence  of  II,j  on  the  evolution  of  the  Reynolds  stress  tensor 
in  (12)  has  motivated  researchers  to  rely  on  homogeneous  turbulent  flows  for  the  testing  and 
calibration  of  pressure-strain  models. 

The  fluctuating  pressure  p  is  a  solution  of  the  Poisson  equation 


^duj  dv, 
dx,  dxj 


du,  duj 
dxj  dx. 
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dui  duj 
dxj  d x , 


(14) 
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which  is  obtained  by  subtracting  the  divergence  of  (4)  from  the  divergence  of  (1).  In  the 
absence  of  boundaries,  equation  (14)  has  the  general  solution 

47 r  7-oo  lx  —  x'|  I  dx',  dx'-  dx',  dx'.  dx'-  dx',  )  ' 


For  homogeneous  turbulent  flows  (where  the  mean  velocity  gradients  are  spatially  uni¬ 
form)  the  pressure-strain  correlation  takes  the  form 


where 


ntJ  —  Aij  -(-  Miju- 


A  1  f°°  ( dux  d2u'ku't  dV' 

,J  47 r  J- oo  ^  dxt )  dx'kdx\  |x'  —  x| 

1  f°°  ( dux  du't  dV' 

,lke  2 tt  7-oo  ^  dxt  )  dx'k  |x'  —  : 


It  has  been  shown  that  Atj  and  Mtjki  are  functionals  -  over  time  and  wavenumber  space  -  of 
the  energy  spectrum  tensor;  see  Weinstock  (1981,  1982)  and  Reynolds  (1987).  In  a  one  point 
closure,  this  dependence  would  suggest  models  for  A+j  and  Mtjki  that  depend  on  the  history 
of  the  Reynolds  stress  tensor  and  dissipation  rate.  The  simplest  such  models  are  algebraic 
in  form: 

Aij  =  eAj(b)  (19) 


Mljki  =  b) 


where 


bij  =  T”  (21) 

TU 

K  =  -\rkk  (22) 

are  the  anisotropy  tensor  and  turbulent  kinetic  energy.  In  (19)  -  (20),  Atj  and  M.xjkt  can 
only  depend  on  rxj  through  bXJ  since  they  are  dimensionless  tensors  that  vanish  in  the  limit  of 
isotropic  turbulence.  Virtually  all  models  for  the  pressure-strain  correlation  that  have  been 
used  in  connection  with  second-order  closure  models  are  of  the  general  form  of  (19)  -  (20). 
The  use  of  this  hierarchy  of  models  for  general  inhomogeneous  turbulent  flows  is  based  on 
the  assumption  of  a  local  homogeneous  structure. 

The  mean  velocity  gradient  tensor  can  be  decomposed  into  symmetric  and  antisymmetric 
parts  as  follows: 

—  Si  j  +  Uij  (23) 
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where 


?  1  /  drij  dv} 

,J  2  V^Xj  <9x, 

_  1  /  dvt  dvj 

2  \5xj  3x, 


are  the  mean  rate  of  strain  tensor  and  mean  vorticity  tensor,  respectively.  Then,  the  model 
for  the  pressure-strain  correlation  specified  by  (19)  -  (20)  can  be  written  in  the  equivalent 
form 

ny=£/‘,ll(b,s',ir')  (26) 


where 


Stj  —  Sa,  Wjj  —  Wij 


are  the  dimensionless  mean  strain  rate  and  vorticity  tensor  whereas  denotes  the  part  of 
the  function  fX]  that  is  linear  in  the  mean  velocity  gradients  and  traceless.  Form  invariance 
under  a  change  of  coordinates  requires  that  transform  as 

Qf(b,  S',  uT/)Qr  =  f(QbQT,QS'QT,Q^QT)  (28) 

where  Q  is  the  rotation  tensor  (and  QT  is  its  transpose)  which  characterizes  a  change  in 
orientation  of  the  coordinate  axes.  In  mathematical  terms,  (28)  requires  that  fx,  be  an 
isotropic  tensor  function  of  its  arguments.  By  using  known  representation  theorems  for 
isotropic  tensor  functions  (see  Smith  1971)  to  construct  fij  -  and  then  by  taking  the  linear 
and  traceless  part  of  ftj  -  the  following  model  for  II, 7  is  obtained: 

n„  =  /?i  ebij  +  02  £(bikbkj  —  -6mn6mn<ftJ)  +  03K~St, 

2 

+P^K(bikSjk  d"  bjlc^ik  ~  X bmn^ rnn^xj) 

3  (29) 

_  _  2  — 

+05 K(bikbiaSjt  +  bjkbiaStt  -  g  b(m  bmn  Srubij  ) 


+P6K(bikWjk  +  bjkUxk)  +  07K(blkbklu>jt  +  bjkbiuu>ti) 


where 


0i  =  Ao(J/,  HI)  +  Ai  (//,  1 1 1)—tr(b  ■  S)  +  0i3{1 1, 1 1 1)—tr(b2  •  S),  i  =  1,2  (30) 

e  e 

^  =  0,(11,111),  j  =  3, 4, ....  7  (31) 

II  =  btjbl},  III  =  btjbjkbkt  (32) 
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arid  tr(-)  denotes  the  trace  (see  Appendix  A  for  more  details).  Equation  (29)  represents  the 
most  general  form  of  the  hierarchy  of  models  (19)  -  (20)  for  the  pressure- strain  correlation 
that  is  consistent  with  the  crucial  physical  constraint  of  invariance  under  coordinate  trans¬ 
formations.  It  will  be  shown  later  that  the  Launder,  Reece,  and  Rodi  model  is  the  linear 
limit  of  (29)  wherein  fix,fiz,fi\,  and  fi$  are  constants  while  /32 , /3s ,  and  fir  are  zero. 

Finally,  in  regard  to  the  general  model,  a  few  comments  should  be  made  concerning  non- 
inertial  frames  of  reference.  In  a  non-inertial  frame,  the  mean  vorticity  tensor  must  be 
replaced  by  the  intrinsic  (i.e.,  absolute)  mean  vorticity  tensor  defined  by 

Wij  =  ZJij  -(-  emjinm  (33) 


where  Clm  is  the  rotation  rate  of  the  non-inertial  frame  relative  to  an  inertial  framing  and  emji 
is  the  permutation  tensor  (see  Launder,  Tselepidakis,  and  Younis  1987  and  Speziale  1989a). 
Furthermore,  Coriolis  terms  must  be  added  to  the  Reynolds  stress  transport  equation  which 
then  takes  the  form 


'  2 

fij  ~  ^ik  Tjk  "h  2  (Tjfcemfcjflm  4-  Tjfcemfcjfim)  (34) 


in  an  arbitrary  non-inertial  reference  frame. 


3.  PLANE  HOMOGENEOUS  TURBULENCE 

We  will  consider  the  general  class  of  plane  homogeneous  turbulent  flows  for  which  the 
mean  velocity  gradient  tensor  can  be  written  in  the  form 

S  -f  u>  0  \ 

0  0  .  (35) 

0  0  / 

Since  the  mean  continuity  equation  (5)  requires  that  dvi/dxi  =  — <5u2/Sx2  in  plane  homoge¬ 
neous  turbulent  flows,  equation  (35)  results  by  simply  aligning  the  coordinates  at  a  45°  angle 
relative  to  the  principal  directions  of  the  symmetric  part  of  dvi/dxj.  Of  course,  in  order  to 
maintain  the  homogeneity  and  two-dimensionality  of  the  mean  flow,  S  and  u  are  constants 
while  Ui  is  given  by 

^  =  (0,0,0)  (36) 

(hence,  the  rotation  is  in  the  plane  where  the  mean  velocity  gradients  are  applied).  Equations 
(35)  -  (36)  encompass,  as  special  cases,  plane  shear,  plane  strain,  rotating  plane  shear,  and 
rotating  plane  strain  turbulence. 


dvi 

dxj 


I  0 

S  —  u 
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The  Reynolds  stress  transport  equation  (34)  for  plane  homogeneous  turbulence  can  be 
written  in  terms  of  the  anisotropy  tensor  (given  that  t*  =  St  is  the  dimensionless 

time)  as  follows: 

§  =  M 


+2—e3kj)  -  bjk{Stk  +  v*k  +  2— e3fcj) 


(37) 


where 


/ 
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(38) 

(39) 


n;  =  n„/2KS 

and  V  ~  Tijdvi/dxj  is  the  turbulence  production.  Equation  (37)  must  be  supplemented 
with  a  transport  model  for  the  turbulent  dissipation  rate  in  order  to  obtain  a  closed  system. 
We  will  consider  the  most  standard  form  of  the  modeled  dissipation  rate  transport  equation 
given  by 

dvi 


.  e  OVi  t 
C‘1KT',dx,  0,2  K 


(40) 


where  CeX  and  Ce2  can  be  functions  of  the  invariants  II  and  III  of  the  anisotropy  tensor 
(in  the  most  commonly  used  form  of  this  model,  Ce i  and  Ce2  are  constants  that  assume  the 
values  of  1.44  and  1.92,  respectively;  see  Hanjalic  and  Launder  1972).  It  will  be  argued  later 
that  some  of  the  crucial  conclusions  to  be  drawn  concerning  the  limitations  of  this  hierarchy 
of  closure  models  for  the  pressure- strain  correlation  are  independent  of  the  specific  form  of 
(40).  Furthermore,  it  should  be  noted  that  virtually  all  of  the  alterations  to  (40)  that  have 
been  proposed  during  the  last  decade  are  highly  ill-behaved  (see  Speziale  1989b). 

The  modeled  dissipation  rate  equation  (40)  can  be  written  in  the  dimensionless  form 


*(■&)-(■&)'*<!*  ~1'7 -to-1* 


(41) 


Equation  (41)  is  obtained  by  combining  (40)  with  the  transport  equation  for  the  turbulent 
kinetic  energy 

K  =  V  -  e  (42) 

which  is  exact  for  homogeneous  turbulence.  When  the  general  model  for  the  pressure-strain 
correlation  (given  by  (29),  with  57^-  replaced  by  W ij)  is  substituted  into  (37),  a  closed  system 
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for  the  determination  of  btj  and  e/ SK  is  obtained.  This  system  of  equations  has  equilibrium 
solutions  of  the  form 

. ...  (j,  |)  (.3) 

(w).-(I'l)  «*> 

where  (•),*,  denotes  the  solution  in  the  limit  as  time  t  — »  oo;  these  solutions  were  shown  by 
Speziale  and  Mac  Giolla  Mhuiris  (1989a)  to  attract  all  initial  conditions.  The  equilibrium 
solutions  (43)  -  (44)  are  obtained  by  solving  the  nonlinear  algebraic  equations  that  result 
when  the  time  derivatives  on  the  left-hand-sides  of  (37)  and  (41)  are  set  to  zero.  It  is  a 
simple  matter  to  show  that  there  is  a  trivial  equilibrium  solution  where 

(s*L=° 

which  exists  for  all  u/S  and  fl/S.  Non-trivial  equilibrium  solutions  where  ( e/SK )0O  ^  0 
exist  for  intermediate  ranges  of  u>jS  and  Q/S  wherein  the  trivial  equilibrium  solution  (45) 
typically  becomes  unstable. 

We  will  now  show  that  the  non-trivial  equilibrium  values  of  ///oo,  (633)00,  and 
(*7«)  OO  are  universal  (i.e.,  completely  independent  of  w/S  and  f l/S)  for  this  hierarchy  of 
models  in  plane  homogeneous  turbulent  flows.  A  closed  system  of  equations  for  the  deter¬ 
mination  of  the  temporal  evolution  of  II,  III,  633  and  e/SK  are  as  follows: 

dll 
df 


dill 
dt * 


2e 

SK 


(>-7)"- 


e  2  V  e 

4033  — 7T77  +  ~‘ 


e SK  3  e SK 


+As*//  +  fts*w-5A7s7? 


(46) 


„  ,  V  e  \  n  TTV  e 
sk  ~  2^U1~SK 


3e 

SK 


(l-j^jlll  +  \lI--±7  -  IIV  S 


e  SK 


£  SK 


Jk^JkU^Jk 

3  V  £  1  TTV  £ 

+4Ai,“7  sic  “  4A//7  sk 


(47) 


1  V  £  1  V  £ 

+  .ft//6337_  _  WH 
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=  -1-  fi  _  Z)h33  _ 

SiC  V  e )  33  3  £  Stf 

(48) 

+  2/?1SKr633  +  2 ^JK  (&33  ~  377) 

di*  ( 5aO  =  (stf)  KC«i  “  ^7  "  (C«2  “  !)]■  (49) 

Equations  (46)  -  (47)  are  obtained  by  multiplying  (37)  with  b  and  b2,  respectively,  and  then 
taking  the  trace  after  the  model  (29)  for  is  substituted.  The  non-trivial  equilibrium 
solutions  are  then  obtained  from  the  nonlinear  algebraic  equations 


4  -  (7)>'~  +  (l)„ 

-(633)*  (|7  +  +  iftj/j, 

+\m ’33)~  (fL-W'-dL 

+  jft((»(i>33)«.  (7)  -  2*7//“  (l)„  =  0 


-(?), 


+5a 


1 


(^33/ao  ~ 


=  0 


(51) 


(Cel-l)f-)  -  (Ce2  -  1)  =  0  (53) 

which  are  derived  by  setting  the  time  derivatives  on  the  left-hand-sides  of  (46)  -  (49)  to 
zero  and  dividing  by  e/SK.  The  system  of  equations  (50)  -  (53)  will  have  solutions  for 
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IIoo,  1 1  loo,  1.^33)00  and  ('P/^oo  that  are  completely  independent  of  a >/  S  and  £1/ S  -  and  hence 
universal  -  for  plane  homogeneous  turbulent  flows. 

This  universal  equilibrium  structure  of  11^,  ///«,,  (633)00  and  {V/e)oo  will  now  be  utilized 
to  obtain  the  simplest  generic  form  of  (29)  that  has  the  same  equilibrium  structure  as  the 
general  model  in  the  phase  space  of  plane  homogeneous  turbulent  flows.  Due  to  these  four 
universal  invariants,  the  quadratic  terms  in  the  rapid  pressure  strain  correlation  are  not 
linearly  independent  for  plane  homogeneous  turbulent  flows.  This  quadratic  part  11^  of  (29) 
is  defined  as  follows: 

nip  =  &K(bikbklS}t  +  b}kbktSu  -  +  prK{bikbklWjt  +  b3kbkiwu).  (54) 

For  plane  homogeneous  turbulent  flows,  a  straightforward,  although  somewhat  tedious,  cal¬ 
culation  yields  the  relationships 

—  _  2  —  _  _  2  —  ?  ITT— 

blkbklSji  +  bjkbktSu  -  -biublmSmk6ij  =  -b33(bikS]k  +  bjkStk  -  -bklSkl8x})  -  -- — Sxj  (55) 

O  OO  633 

btkbiaW ji  +  bjkbkiWii  —  —b33(btkW jk  -f  b]k W lk )  (56) 

where  we  have  made  use  of  (38)  and  the  fact  that  the  anistropy  tensor  is  of  the  form 


6tJ  = 


611  612  0 

612  622  0 

0  0  633 


(57) 


in  such  flows.  Due  to  (55)  -  (56),  and  the  fact  that  1 /<»,  1 1 1^,  (633)00  and  ( V/e)^  are 
universal  invariants  for  all  plane  homogeneous  turbulent  flows,  ic  follows  that  the  quadratic 
terms  in  the  rapid  pressure-strain  are  directly  related  to  the  linear  terms  in  such  flows. 
Consequently,  the  equilibrium  structure  of  (29)  in  plane  homogeneous  turbulent  flows  will 
be  indistinguishable  from  that  of  the  substantially  simpler  model 


n,j  —  Ciebij  +  c2e(bikbkj  jbmnbmn8jj) 
+c3k~Sij  +  c^K(bikSjk  +  bjk^ik 

mn^mn  Sij)  +  c5K(btkW]k  f  ojkWtk) 


(58) 


-lb, 


where  Ci,c2, ... ,  c5  are  dimensionless  constants  and  we  have  made  use  of  the  fact  that 

Jr(b-5)=-~  (59) 

tr(b2  -5)  =  (60) 
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which  was  also  used  in  the  derivation  of  (46)  -  (49).  In  alternative  terms,  (58)  is  topolog¬ 
ically  equivalent  to  the  general  model  (29)  in  so  far  as  the  equilibrium  structure  of  plane 
homogeneous  turbulent  flows  is  concerned. 

It  is  rather  striking  that  an  analysis  of  the  equilibrium  states  of  arbitrary  plane  ho¬ 
mogeneous  turbulence  -  coupled  with  the  crucial  physical  constraint  of  invariance  under 
coordinate  transformations  -  collapses  the  general  pressure-strain  model 

n a  =  tX3{ b)  +  KMiilu( b)^,  (6i) 

(which  can  have  as  many  as  forty-five  independent  functions  of  b)  to  the  substantially  sim¬ 
plified  model  (58)  that  has  only  five  undetermined  constants.  In  the  next  section,  a  new 
model  for  the  pressure-strain  correlation  will  be  developed. 

4.  THE  SSG  MODEL:  ITS  ASYMPTOTIC  ANALYSIS  AND  CALIBRATION 

Now,  a  new  model  for  the  pressure-strain  correlation  -  which  we  will  call  the  SSG  model 
-  will  be  developed  based  on  the  previous  invariant  dynamical  systems  analysis  coupled  with 
the  following  additional  constraints: 

(i)  Asymptotic  consistency  in  the  limit  of  small  anisotropies 

(ii)  Consistency  with  Rapid  Distortion  Theory  for  homogeneously  strained  turbulent  flows 
that  are  initially  isotropic 

(iii)  Consistency  with  the  equilibrium  values  for  homogeneous  shear  flow  obtained  from  the 
physical  experiments  of  Tavoularis  and  Corrsin  (1981) 

(iv)  Consistency  with  the  RDT  results  of  Bertoglio  (1982),  for  rotating  shear  flows,  which 
predict  that  the  most  unstable  flow  occurs  when  the  ratio  of  the  rotation  rate  to  the 
shear  rate  0/ S  =  0.25  and  that  a  flow  restabilization  occurs  when  f 1/ S  >  0.5 

(v)  Consistency  with  the  results  of  physical  experiments  on  the  decay  of  isotropic  turbu¬ 
lence  and  the  return  to  isotropy  of  an  initially  anisotropic,  homogeneous  turbulence. 

Since  the  magnitude  of  the  anisotropy  tensor  is  relatively  small  (||b||  <  0.20  for  most  turbu¬ 
lent  flows  of  engineering  and  scientific  interest),  we  feel  that  terms  which  are  of  a  comparable 
order  of  magnitude  in  bxj  should  be  maintained  unless  there  is  some  overriding  physical  rea¬ 
son  not  to  do  so.  In  this  fashion,  the  model  can  then  be  thought  of  as  an  asymptotically 
consistent  truncation  of  a  Taylor  series  expansion  of  AXJ  and  in  the  variable  biy  Since 
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the  simplified  model  for  the  rapid  pressure-strain  correlation  in  (58)  is  of  0(b),  this  suggests 
that  c3  -  which  in  general  can  be  a  function  of  the  invariants  of  b  -  should  be  replaced  by 

c3  -  c;ii* 

(where  C3  and  are  constants)  for  asymptotic  consistency.  Furthermore,  since  the  model 
for  the  slow  pressure-strain  correlation  is  of  0(b2),  and  since  most  engineering  flows  have 
significant  regions  where  V  >  e,  we  will  replace  the  constant  Cj  with  the  coefficient 


-(Ox  +  C[V/e) 


where  C\  and  C*  are  constants.  This  yields  the  following  model  for  the  pressure-strain 
correlation: 


IT*. j  —  — (C  is  +  C*V)bij  +  CzEibikbkj 


-\bmnbmn6X])  +  (C3  -  Cllti)KSij 

_  _  2  — 
+C<K(btkSjk  +  b]kStk  -  ~b  rnn  S mn^ij  ) 


(62) 


+CsK(blkW,k  +  4*^*). 

Although  (62)  is  topologically  equivalent  to  (58)  in  so  far  as  the  equilibrium  states  are 
concerned,  it  will  give  rise  to  different  temporal  evolutions.  We  feel  that  it  is  better  to 
use  (62)  as  our  final  model  for  the  pressure-strain  correlation  since  all  terms  that  are  of  a 
comparable  order  in  6tJ  have  been  maintained  for  asymptotic  consistency. 

Before  using  constraints  (ii)  -  (v)  to  calibrate  the  SSG  model  given  by  (62),  a  few  com¬ 
ments  are  in  order  concerning  the  relationship  between  this  new  model  and  previously  pro¬ 
posed  models.  The  SSG  model  is  not  significantly  more  complicated  than  the  commonly 
used  Launder,  Reece,  and  Rodi  model  which  can  be  obtained  from  (62)  in  the  linear  limit  as 
C*,C2  and  CZ  go  to  zero.  In  fact,  the  SSG  model  is  substantially  simpler  than  the  recently 
proposed  nonlinear  models  of  Shih  and  Lumley  (1985)  and  Fu,  Launder,  and  Tselepidakis 
(1987)  (see  Appendix  B). 

The  coefficients  C\  and  C3  have  been  calibrated  by  considering  the  return  to  isotropy 
problem  (see  Sarkar  and  Speziale  1989).  Of  course,  for  the  return  to  isotropy  problem,  only 
the  terms  containing  the  coefficients  C\  and  C2  in  the  pressure- strain  correlation  survive  (i.e., 
the  rapid  pressure-strain  correlation  vanishes).  Based  upon  realizability,  dynamical  systems 
considerations,  and  the  phase  space  portrait  of  return  to  isotropy  experiments,  the  following 
values  of  C\  and  C3  were  arrived  at  by  Sarkar  and  Speziale  (1989): 


Cx  =  3.4 


(63) 
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C2  =  3(Cj  -  2)  =  4.2.  (64) 

Interestingly  enough,  the  value  of  C\  =  3.4  is  quite  close  to  the  value  of  3.6  for  the  Rotta 
coefficient  that  is  currently  used  in  the  basic  model  of  Launder  and  his  co-workers.  However, 
as  demonstrated  by  Sarkar  and  Speziale  (1989),  the  quadratic  term  containing  C2  is  crucial 
to  properly  capture  the  experimental  trends.  In  Figure  1,  the  predictions  of  the  SSG  model 
and  the  Launder,  Reece,  and  Rodi  (LRR)  model  are  compared  in  the  £  —  rj  phase  space 
with  the  experimental  data  of  Choi  and  Lumley  (1984)  for  the  return  to  isotropy  from  plane 
strain.  The  SSG  model  exhibits  a  curved  trajectory  that  is  well  within  the  range  of  the 
experimental  data;  the  LRR  model  -  as  well  as  any  model  for  which  C2  =  0  -  erroneously 
predicts  a  straight  line  trajectory.  In  Figures  2(a)  -  (b),  the  predictions  of  the  SSG  model 
and  the  LRR  model  for  the  temporal  evolution  of  the  anisotropy  tensor  are  compared  with 
experimental  data  for  the  relaxation  from  plane  strain  experiment  of  Choi  and  Lumley  (1984) 
and  plane  contraction  experiment  of  Le  Penven,  Gence,  and  Comte-Bellot  (1985).  The  SSG 
model,  on  balance,  yields  improved  predictions  over  the  LRR  model.  For  more  detailed 
discussions  and  comparisons,  the  reader  is  referred  to  the  paper  by  Sarkar  and  Speziale 
(1989)  where  this  quadratic  model  for  the  slow  pressure-strain  correlation  was  compared 
with  data  from  four  independent  experiments  on  the  return  to  isotropy. 

Constraint  (ii),  which  requires  consistency  with  RDT  for  a  homogeneously  strained  tur¬ 
bulence  that  is  initially  isotropic,  is  commonly  enforced  in  the  turbulence  modeling  literature. 
While  the  dynamical  systems  analysis  presented  in  Section  3  can  guarantee  proper  long-time 
behavior,  this  RDT  constraint  can  be  of  considerable  assistance  in  ensuring  proper  short- 
time  behavior;  if  a  model  properly  captures  both  the  short  and  long-time  behavior,  it  stands 
an  excellent  chance  of  performing  well  for  all  times.  This  RDT  result  requires  that  (see  Crow 
1968) 

lim  n„  =  \k~Sij  (65) 

b— >o  5 

and,  hence,  that 

C,  =  j  (66) 

for  the  SSG  model.  We  found  that  models  which  deviated  significantly  from  (66)  performed 
poorly  in  homogeneous  shear  flows  (e.g.,  such  models  yielded  spurious  points  of  inflection  in 
the  time  evolution  of  the  turbulent  kinetic  energy). 

Constraints  (iii)  -  (iv)  were  used  to  calibrate  the  remaining  constants  in  the  model, 
namely,  C*,  C3,  C4,  and  C5  as  well  as  the  constant  Cel  in  the  modeled  e-transport  equation. 
This  was  done  using  a  value  of  Ce2  —  1.83  (as  opposed  to  the  more  commonly  adopted 
value  of  1.92)  which  yields  a  power  law  decay  in  isotropic  turbulence  with  an  exponent  of 
1.2  -  a  value  which  is  in  better  agreement  with  available  experimental  data  as  discussed  by 


13 


Reynolds  (1987).  It  was  not  possible  to  obtain  the  exact  equilibrium  values  of  Tavoularis 
and  Corrsin  (1981)  for  homogeneous  shear  flow  and  satisfy  the  RDT  results  of  Bertoglio  for 
rotating  shear  flow  (i.e.,  constraints  (iii)  -  (iv)).  A  numerical  optimization  yielded  the  values 
of  C *  =  1.80,  =  1.30,  C4  =  1.25,  C5  =  0.40  and  Ce  1  =  1.44  as  the  best  compromise.  The 

equilibrium  values  obtained  from  the  SSG  model  (using  these  values  of  the  constants)  for 
homogeneous  shear  flow  are  compared  with  the  values  obtained  from  the  Launder,  Reece, 
and  Rodi  model  and  the  experiments  of  Tavoularis  and  Corrsin  (1981)  in  Table  1.  From 
these  results,  it  is  clear  that  the  predictions  of  the  SSG  model  are  well  within  the  range 
of  the  experiments  whereas  the  predictions  of  the  Launder,  Reece,  and  Rodi  model  deviate 
significantly.  Furthermore,  the  SSG  model  predicts  that  the  largest  growth  rate  in  rotating 
shear  flow  occurs  when  0/5  ~  0.22  and  that  a  flow  restabilization  occurs  when  0/5  >  0.53 
in  comparison  to  the  corresponding  RDT  results  of  0/5  =  0.25  and  0/5  >  0.5  (see  Bertoglio 
1982).  These  predictions  of  the  SSG  model  are  considerably  better  than  those  of  the  LRR 
model  which  erroneously  predicts  that  the  largest  growth  rate  occurs  when  0/5  =  0.14 
and  that  a  flow  restabilization  occurs  when  0/5  >  0.37.  A  more  detailed  discussion  of  the 
performance  of  the  models  in  rotating  shear  flow  will  be  presented  in  the  next  section. 

The  SSG  model  has  been  carefully  calibrated  to  perform  well  in  shear  flows  both  with 
and  without  added  rotational  strains.  It  is  our  belief  that  this  will  significantly  enhance 
the  performance  of  the  model  in  turbulent  boundary  layers  with  streamline  curvature  -  an 
analogous  flow  with  a  variety  of  important  engineering  applications.  However,  unlike  other 
recently  derived  models  for  the  pressure-strain  correlation  such  as  the  Shih-Lumley  (1985) 
model  and  the  Fu,  Launder,  and  Tselepidakis  (1987)  model,  the  SSG  model  given  by  equation 
(62)  does  not  satisfy  the  strong  form  of  realizability.  The  strong  form  of  realizability  (see 
Lumley  1978)  constitutes  a  sufficient  condition  to  guarantee  positive  component  energies 
in  homogeneous  turbulent  flows.  The  SSG  model  only  satisfies  a  weak  form  of  realizability 
wherein  the  turbulent  kinetic  energy  is  guaranteed  to  be  positive;  this  is  a  direct  consequence 
of  the  form  of  the  modeled  f-transport  equation  (see  Speziale  1989b).  We  decided  to  opt  for 
the  weaker  form  of  realizability  for  two  major  reasons.  First,  if  the  turbulent  kinetic  energy 
is  positive,  realizability  can  only  be  violated  by  fairly  large  anisotropies,  such  that 

IN  >  i 

(where  ||  •  ||  is  the  L2  norm  or  maximum  eigenvalue)  which  are  outside  of  the  expected  domain 
of  applicability  of  such  idealized  models.  It  must  be  kept  in  mind  that,  so  long  as  the  model 
yields  a  positive  turbulent  kinetic  energy,  it  can  be  applied  to  a  flow  (it  is  primarily  negative 
kinetic  energies  that  are  computationally  fatal).  Second,  it  has  been  our  experience  that 
models  which  satisfy  the  strong  form  of  realizability  become  computationally  “stiff”  in  flows 
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with  large  anisotropies.  This  results  from  the  fact  that  the  finite  difference  form  of  the 
modeled  equations  usually  do  not  exactly  satisfy  realizability  (see  Speziale  and  Mac  Giolla 
Mhuiris  1989a).  No  such  problems  are  encountered  by  the  weak  form  of  realizability  since  it  is 
satisfied  exactly  by  most  standard  numerical  formulations  of  the  model.  Finally,  it  should  be 
mentioned  that  the  SSG  model  was  not  forced  to  satisfy  material  frame  indifference  (MFI)  in 
the  limit  of  two-dimensional  turbulence  (Speziale  1983)  which  constitutes  another  extreme 
constraint  that  is  a  rigorous  consequence  of  the  Navier-Stokes  equations.  It  has  recently 
become  apparent  to  us  that  when  such  extreme  constraints  as  MFI  and  strong  realizability 
(correct  as  they  may  be  for  general  flows)  are  applied  to  highly  idealized  models,  there  is  a 
strong  possibility  that  the  model  will  become  overly  biased  so  that  it  performs  poorly  in  the 
more  commonly  encountered  turbulent  flows. 


5.  PERFORMANCE  OF  THE  SSG  MODEL  IN  HOMOGENEOUS  FLOWS 


The  SSG  model  given  by  equation  (62)  will  now  be  tested  in  four  independent  homo¬ 
geneous  turbulent  flows.  For  the  purpose  of  clarity,  we  will  summarize  the  values  of  the 
constants  that  were  arrived  at  in  the  previous  section: 


cx  =  3.4,  C{  =  1.80,  C2  =  4.2 

(67) 

Co  =  7,  C;  =  1.30,  C4  =  1.25 

5 

(68) 

5  =  0.40,  Cei  =  1-44,  Ce2  =  1-83. 

(69) 

The  problem  of  homogeneous  turbulent  shear  flow  in  a  rotating  frame  will  be  considered 
first.  For  this  case,  the  mean  velocity  gradients  and  the  rotation  rate  of  the  reference  frame 
are  given  by 


dvl 

dxj 


(  0  S  0  \ 
0  0  0 
V  o  o  o  / 


n.  =  (o.o.fi), 


(70) 

(71) 


respectively,  in  matrix  form.  The  initial  conditions  correspond  to  a  state  of  isotropic  turbu¬ 
lence  where 

bij  =  0,  K  =  Kq ,  e  —  Co  (72) 

at  time  t  =  0.  It  was  shown  by  Speziale  and  Mac  Giolla  Mhuiris  (1989a)  that  the  solution 
only  depends  on  the  initial  conditions  through  the  dimensionless  parameter  Co/SKq]  the 
dependence  of  the  solution  on  the  shear  rate  and  rotation  rate  is  exclusively  through  the 
dimensionless  parameter  0,/S.  Two  types  of  equilibrium  solutions  have  been  established 
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for  this  problem  (Speziale  and  Mac  Giolla  Mhuiris  1989a):  one  where  ( e/SK )ao  =  0  which 
exists  for  all  0/5  and  one  where  ( e/SK)00  >  0  which  exists  only  for  a  small  intermediate 
band  of  values  for  f 2/5.  The  zero  equilibrium  value  is  associated  predominantly  with  stable 
flow  wherein  K  and  £  undergo  a  power-law  time  decay;  the  nonzero  equilibrium  values  are 
associated  with  unstable  flow  wherein  K  and  £  undergo  an  exponential  time  growth.  The  two 
solutions  undergo  an  exchange  of  stabilities  for  intermediate  values  of  f 1/5  (which  includes 
the  case  of  pure  shear  where  0/5  =  0).  In  this  fashion,  the  second-order  closures  are  able 
to  account  for  both  the  shear  instability  -  with  its  exponential  time  growth  of  disturbance 
kinetic  energy  -  and  the  stabilizing  (or  destabilizing)  effect  of  rotations  on  shear  flow. 

In  Figures  3(a)  -  (c),  the  predictions  of  the  SSG  model  for  the  time  evolution  of  turbulent 
kinetic  energy  are  compared  with  those  of  the  LRR  model  and  the  results  of  the  large-eddy 
simulations  of  Bardina,  Ferziger,  and  Reynolds  (1983)  for  rotating  shear  flow.  From  these 
figures,  it  is  clear  that  the  SSG  model  does  a  much  better  overall  job  of  capturing  the  trends 
of  the  large-eddy  simulations.  Several  observations  are  noteworthy:  (a)  the  LRR  model 
exhibits  too  strong  a  growth  rate  for  pure  shear  (f l/S  =  0)  in  comparison  to  the  SSG  model 
and  large-eddy  simulations,  (b)  both  the  SSG  and  LRR  models  yield  too  weak  a  growth 
rate  for  the  0/5  =  0.25  case,  however,  the  SSG  model  is  substantially  better,  and  (c)  the 
SSG  model  properly  captures  the  weak  growth  rate  that  occurs  for  0/5  =  0.5,  whereas  the 
LRR  model  erroneously  predicts  a  flow  restabilization.  The  premature  flow  restabilization 
predicted  by  the  LRR  model  at  0/5  %  0.37  is  somewhat  serious  since,  in  addition  to  the 
results  of  large-eddy  simulations,  linear  stability  theory  and  RDT  predict  that  there  should 
be  unstable  flow  for  the  entire  range  of  0  <0/5  <  0.5  (see  Lezius  and  Johnston  1976  and 
Bertoglio  1982).  As  mentioned  earlier,  the  SSG  model  does  not  predict  a  flow  restabilization 
until  0/5  >  0.53. 

It  would  be  useful  at  this  point  to  compare  the  performance  of  the  SSG  model  in  rotating 
shear  flow  with  that  of  some  newer  models  that  have  been  recently  proposed.  Three  such 
models  -  the  model  of  Shih  and  Lumley  (1985),  the  model  of  Fu,  Launder,  and  Tselepidakis 
(1987),  and  the  RNG  model  of  Yakhot  and  Orszag  (1988)  -  were  compared  in  a  recent  study 
of  Speziale,  Gatski,  and  Mac  Giolla  Mhuiris  (1989).  It  was  established  in  that  study  that 
the  Fu,  Launder,  and  Tselepidakis  (FLT)  model  performed  the  best  among  these  models 
in  rotating  shear  flow.  Hence,  for  simplicity,  we  will  only  compare  the  SSG  model  with 
the  FLT  model.  In  Figures  4(a)  -  (c),  the  results  for  the  turbulent  kinetic  energy  obtained 
from  the  SSG  model  and  the  FLT  model  for  the  rotation  rates  of  0/5  =  0,0/5  =  0.25, 
and  0/5  =  0.50  are  compared  with  the  large-eddy  simulations  of  Bardina,  Ferziger,  and 
Reynolds  (1983)  for  rotating  shear  flow.  It  is  clear  from  these  results  that  the  SSG  model 
properly  captures  the  trends  of  the  large-eddy  simulations  which  indicate  that  all  three  cases 
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are  unstable  and  that  the  Cl/S  =  0.25  case  has  the  strongest  growth  rate.  On  the  otherhand, 
the  FLT  model  erroneously  predicts  that  the  Cl/S  =  0  and  Cl/S  =  0.25  cases  are  equally 
energetic  and  that  the  Cl/S  =  0.5  case  has  undergone  a  restabilization.  Comparable  to  the 
LRR  model,  the  FLT  model  erroneously  predicts  a  premature  restabilization  at  Cl /  S  ~  0.39. 
It  may  be  of  concern  that  a  heavy  emphasis  has  been  placed  on  comparisons  with  large- 
eddy  simulations  for  rotating  shear  flow  (unfortunately,  no  direct  simulations  or  physical 
experiments  have  been  conducted  for  this  problem).  However,  it  must  be  emphasized  that 
the  critical  evaluations  have  been  based  on  which  states  should  be  more  energetic  -  results 
which  have  been  confirmed  independently  by  RDT  and  linear  stability  theory. 

A  bifurcation  diagram  for  the  general  hierarchy  of  closure  models  (61)  is  shown  in  Figure 
5  for  rotating  shear  flow.  Here,  the  equilibrium  value  of  (e/SK)^  is  plotted  as  a  function  of 
Cl /  S.  The  SSG  model  as  well  as  the  other  commonly  used  models  have  the  same  topological 
structure  in  rotating  shear  flow  as  indicated  in  Figure  5.  There  are  two  equilibrium  solutions: 
the  solution  where  (e/ SK)00  =  0  exists  for  all  Cl/S  but  becomes  unstable  in  the  interval  AB\ 
the  nonzero  solution  for  (e/SK)oa,  which  lies  on  the  semi-ellipse  ACB,  exchanges  stabilities 
with  the  trivial  solution  (e/5Ar)00  =  0  in  the  interval  A  <  Cl/S  <  B.  For  Cl/S  <  A  —  SA  and 
Cl/ S  >  B  A  SB  (where  8 A  and  SB  represent  a  small  increment  that  depends  on  the  model) 
the  trivial  equilibrium  value  of  (e/ SK)00  —  0  is  associated  with  solutions  where  the  kinetic 
energy  undergoes  a  power  law  decay  with  time;  for  A  —  8A  <  Cl/S  <  B  +  SB,  this  trivial 
solution  is  associated  with  solutions  where  the  kinetic  energy  undergoes  a  power  law  growth 
with  time.  The  nonzero  equilibrium  values  (e/SK)00  >  0  (on  the  semi-ellipse  ACB)  are 
associated  with  solutions  where  the  kinetic  energy  grows  exponentially  with  time.  It  can  be 
shown  (see  Speziale  and  Mac  Giolla  Mhuiris  1989a)  that  the  growth  rate  A  for  A  <  Cl/ S  <  B 
is  given  by 

(73) 

where  a  =  (Ce2  —  l)/(Cei  —  1).  Hence,  point  C  -  which  corresponds  to  the  maximum  value 
of  (e/5/£r)o0  -  represents  the  most  energetic  state  with  the  largest  growth  rate  of  kinetic 
energy. 

The  coordinates  [Cl/S,  (e/SK)^}  of  points  A,  B,  and  C  (in  Figure  5)  for  the  Launder, 
Reece,  and  Rodi  model  and  the  SSG  model  are  given  below: 

LRR  Model 

A  =  [-0.09,0],  B  =  [0.37,0],  C  =  [0.14,0.167] 


SSG  Model 


A  =  [-0.09,0],  B  =  [0.53,0],  C  =  [0.22,0.254], 
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The  improved  performance  of  the  SSG  model  in  rotating  shear  flow  is  largely  due  to  the 
fact  that  its  most  energetic  state  (point  C  on  the  bifurcation  diagram  shown  in  Figure  5)  is 
located  in  close  proximity  to  Cl/S  =  0.25  -  the  value  predicted  by  rapid  distortion  theory. 
However,  it  needs  to  be  mentioned  at  this  point  that  the  reason  we  were  not  able  to  satisfy 
this  RDT  result  exactly  is  due  to  a  defect  in  the  general  hierarchy  of  models  (61).  Due  to 
(73)  and  the  fact  that  the  bifurcation  diagram  is  symmetric  about  its  most  energetic  state 
(point  C  in  Figure  5),  the  general  hierarchy  of  models  (61)  erroneously  predicts  Richardson 
number  similarity  if  point  C  is  located  at  Cl/S  =  0.25.  Such  models  will  yield  solutions  for 
K  and  t  that  scale  with  the  Richardson  number 

Ri  =  -2(Q/5)(1  -  20/5)  (74) 

and,  thus,  erroneously  predict  that  the  Cl/S  =  0  and  Cl/S  =  0.5  cases  are  identical.  Large- 
eddy  simulations,  RDT  (Bertoglio  1982),  and  independent  mathematical  analyses  of  the 
Navier-Stokes  equations  (Speziale  and  Mac  Giolla  Mhuiris  1989b)  indicate  that  the  fi/5  =  0 
and  Cl/S  =  0.5  cases  are  distinct.  By  moving  the  most  energetic  state  a  small  distance  to 
the  left  of  Cl/S  =  0.25  -  as  is  done  with  the  SSG  model  -  the  proper  growth  rates  obtained 
from  large-eddy  simulations  for  Cl/S  —  0  and  Cl/S  —  0.5  can  be  reproduced  (see  Figures  3(a) 
and  3(c)).  However,  the  substantially  larger  growth  rate  for  Cl/S  =  0.25  shown  in  Figure 
3(b)  (which  has  independent  support  in  the  RDT  calculations  of  Bertoglio  1982),  cannot  be 
reproduced  by  the  SSG  model.  This  is  a  defect  in  the  general  hierarchy  of  models  (61)  that  is 
intimately  tied  to  their  prediction  of  universal  equilibrium  values  for  //<», ///qo,  (633)^  and 
(Ty£)oo  in  plane  homogeneous  turbulent  flows  -  an  oversimplification  that  is  not  supported 
by  physical  or  numerical  experiments. *  Nonetheless,  despite  this  deficiency,  the  SSG  model 
performs  reasonably  well  -  and  is  superior  to  other  existing  second-order  closures  -  for 
rotating  shear  flow  as  evidenced  by  Table  1,  Figures  3(a)  -  (c),  and  Figures  4(a)  -  (c). 

Now,  we  will  examine  the  performance  of  the  SSG  model  in  homogeneous  plane  strain 
turbulence  for  which  the  mean  velocity  gradients  are  given  by 

o-  (  S  0  0  \ 

0  -5  0  (75) 

dX}  \  0  0  0  / 

and  the  turbulence  evolves  from  an  initial  state  of  isotropy.  Comparisons  of  the  model  predic¬ 
tions  will  be  made  with  the  direct  numerical  simulations  of  Lee  and  Reynolds  (1985)  on  plane 
strain.  Such  comparisons  must  be  made  with  caution  due  to  the  low  turbulence  Reynolds 
numbers  of  the  direct  simulations.  However,  comparisons  with  physical  experiments  (e.g., 


tit  is  not  possible  to  tie  this  deficiency  to  the  modeled  e-transport  equation  since  all  dependence  on  e 
can  be  eliminated  in  the  RDT  limit. 
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Tucker  and  Reynolds  1968  and  Gence  and  Mathieu  1980)  are  equally  problematical  due  to 
the  uncertainty  in  the  initial  conditions  for  c/SK  and  possible  large-scale  contamination 
from  the  walls  of  the  test  apparatus. 

In  Figure  6,  the  time  evolution  of  the  turbulent  kinetic  energy  for  the  LRR  model  and  SSG 
model  are  compared  with  the  direct  simulations  of  Lee  and  Reynolds  (1985)  for  plane  strain 
corresponding  to  the  initial  condition  Cq/SKo  =  2.0.  It  is  clear  from  this  figure  that  both 
models  perform  extremely  well.  In  Figure  7,  the  time  evolution  of  the  non-zero  components 
of  the  anisotropy  tensor  are  shown.  Although  the  quantitative  accuracy  of  the  models  is 
not  extremely  good,  it  is  clear  that  the  SSG  model  does  better  than  the  LRR  model  and 
reproduces  the  crucial  trends  of  the  direct  simulations.  In  Figures  8  -  9,  the  time  evolution 
of  the  turbulent  kinetic  energy  and  non-zero  components  of  the  anisotropy  tensor  are  shown 
corresponding  to  the  initial  condition  of  Co/SKo  =  1.0.  The  same  conclusions  can  be  drawn 
from  these  results:  the  SSG  model  yields  improved  predictions  over  the  LRR  model  and, 
on  balance,  compares  reasonably  well  with  the  direct  simulations  which  would  be  expected 
to  have  somewhat  elevated  anisotropies  due  to  the  lower  turbulence  Reynolds  number.  We 
will  not  make  more  extensive  comparisons  with  the  predictions  of  other  turbulence  models 
since  our  main  purpose  here  was  to  simply  establish  that  the  alterations  made  in  the  LRR 
model  -  to  yield  the  SSG  model  with  its  improved  behavior  in  rotating  shear  flows  -  do  not 
compromise  its  performance  in  plane  strain. 

Finally,  we  will  examine  the  performance  of  the  SSG  model  for  the  axisymmetric  contrac¬ 
tion  and  expansion  in  homogeneous  turbulence.  Since  the  SSG  model  (like  virtually  all  other 
existing  models  for  the  pressure-strain  correlation)  was  calibrated  based  on  plane  homoge¬ 
neous  turbulent  flows,  it  would  be  desirable  to  assess  its  performance  in  a  three-dimensional 
flow.  For  the  axisymmetric  contraction,  the  mean  velocity  gradients  are  given  by 


dvt 


dx} 


(SO  0  \ 

0  -i  S  0 

\0  0  ~\s ) 


whereas  in  the  axisymmetric  expansion  they  take  the  form 


dvi 


dxj 


-S  0  0  \ 

0  iS  0 

0  0  \S  } 


(76) 


(77) 


The  time  evolution  of  each  of  these  turbulent  flows  -  from  an  initially  isotropic  state  -  will 
be  considered.  Hence,  as  with  plane  shear  and  plane  strain,  the  solutions  will  only  depend 
on  the  initial  conditions  through  the  parameter  e0/SK0.  Comparisons  will  be  made  with 
the  predictions  of  the  LRR  model  and  the  direct  numerical  simulations  of  Lee  and  Reynolds 
(1985)  for  the  same  reasons  as  cited  above. 
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In  Figure  10,  the  time  evolution  of  the  turbulent  kinetic  energy  for  the  axisymmetric 
contraction  is  shown  corresponding  to  the  initial  condition  to/ SKq  =  0.179  which  was  taken 
from  the  direct  simulations  of  Lee  and  Reynolds  (1985).  From  these  results,  it  is  clear  that 
the  SSG  model  yields  noticeably  improved  predictions  over  the  LRR  model;  however,  both 
models  predict  growth  rates  that  are  smaller  than  those  in  the  direct  simulations.  In  Figure 
11,  the  time  evolution  of  the  nonzero  components  of  the  anisotropy  tensor  are  shown  for  the 
axisymmetric  contraction  where  e0/SK0  =  0.179.  While  the  differences  between  the  SSG 
and  LRR  models  is  small,  it  is  clear  that  the  SSG  model  yields  results  that  are  more  in  line 
with  the  direct  simulations. 

In  Figure  12,  the  time  evolution  of  the  turbulent  kinetic  energy  for  the  axisymmetric 
expansion  is  shown  for  the  initial  condition  e0/SK0  =  2.45  taken  from  the  direct  numerical 
simulations  of  Lee  and  Reynolds  (1985).  It  is  clear  from  this  figure  that  both  the  SSG 
and  LRR  models  yield  results  that  are  in  excellent  agreement  with  the  direct  simulations. 
However,  the  time  evolution  of  the  nonzero  components  of  the  anisotropy  tensor  shown  in 
Figure  13  show  more  significant  differences.  Here,  the  predictions  of  the  SSG  model  appear 
to  be  substantially  better  than  those  of  the  LRR  model. 

6.  CONCLUSIONS 

The  modeling  of  the  pressure-strain  correlation  of  turbulence  has  been  considered  based 
on  invariance  arguments  and  a  dynamical  systems  approach.  Several  important  conclusions 
have  been  drawn  about  the  standard  hierarchy  of  closures  (61)  which  led  to  the  development 
of  a  new  model  -  the  SSG  model.  A  summary  of  these  findings  can  be  given  as  follows: 

(i)  It  was  proven  that  the  standard  hierarchy  of  models  yields  non-trivial  values  for  the 
equilibrium  states  //<„,  IIIoo,  (633)00,  and  ('P/£)00  that  are  universal  (i.e. ,  that  do  not 
depend  on  cu/5,  0/5  or  the  initial  conditions)  for  plane  homogeneous  turbulent  flows. 
As  a  direct  consequence  of  these  universal  invariants,  it  was  shown  that,  for  plane 
homogeneous  turbulent  flows,  the  general  model  (61)  for  the  pressure-strain  correlation 
is  topologically  equivalent  to  a  substantially  simpler  model  -  the  SSG  model  -  which 
is  only  quadratically  nonlinear  in  the  anisotropy  tensor. 

(ii)  The  SSG  model  was  calibrated  by  using  existing  data  from  isotropic  decay  experiments, 
return  to  isotropy  experiments,  and  homogeneous  shear  flow  experiments  along  with 
the  RDT  results  of  Crow  (1968)  and  Bertoglio  (1982).  By  means  of  this  more  systematic 
method  of  calibration,  the  SSG  model  was  demonstrated  to  perform  better  than  the 
Launder,  Reece,  and  Rodi  model  -  as  well  as  the  newer  models  of  Shih  and  Lumley 
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and  Fu,  Launder,  and  Tselepidakis  -  for  a  variety  of  homogeneous  turbulent  flows. 
The  flows  that  were  examined  included  the  challenging  test  case  of  rotating  shear  flow 
(where  rotations  can  either  stabilize  or  destabilize  the  flow)  and  the  axisymmetric 
expansion/contraction  which  constitutes  a  three-dimensional  mean  turbulent  flow. 

(iii)  Although  the  SSG  model  performs  reasonably  well  for  a  variety  of  homogeneous  tur¬ 
bulent  flows,  there  are  still  major  deficiencies  with  it  that  are  intrinsic  to  this  general 
hierarchy  of  models.  These  deficiencies,  emanate  from  the  prediction  of  universal  equi¬ 
librium  values  for  7/qo,  1 1  loo,  (633)00  and  ( V/e )«,  in  plane  homogeneous  turbulent  flows 
-  an  obvious  oversimplification  that  is  not  supported  by  physical  experiments.  As  a 
result  of  this  deficiency,  the  general  model  (61)  erroneously  predicts  that  rotating  shear 
flow  has  growth  rates  that  are  symmetric  about  their  most  energetic  value.  Hence,  in 
order  to  satisfy  the  RDT  constraint  of  Bertoglio  (1982)  -  which  puts  the  most  energetic 
state  at  Q/S  =  0.25  -  the  models  must  exhibit  Richardson  number  similarity.  This  is 
inconsistent  with  the  Navier-Stokes  equations  as  proven  by  Speziale  and  Mac  Giolla 
Mhuiris  (1989b)  and  illustrated  by  Bardina,  Ferziger,  and  Reynolds  (1983). 

(iv)  Since  the  general  model  (61)  for  the  pressure-strain  correlation  gives  an  incomplete 

picture  of  plane  homogeneous  turbulent  flows  no  matter  what  form  is  taken  for  Aj(b) 
and  b),  we  feel  that  the  process  of  adding  highly  complex  nonlinear  terms  in 

is  somewhat  questionable.  Such  complex  nonlinear  terms  in  the  rapid  pressure-strain 
correlation  have  been  largely  motivated  by  the  desire  to  satisfy  the  strong  form  of 
realizability.  However,  it  must  be  remembered  that  the  strong  form  of  this  constraint 
only  constitutes  a  sufficient  condition  for  the  satisfaction  of  realizability  in  homoge¬ 
neous  turbulent  flows.  Due  to  the  relatively  simple  topological  structure  of  the  general 
model  in  rotating  shear  flow  -  which  is  in  no  way  altered  by  the  addition  of  more 
complex  nonlinearities  in  b^j  -  the  application  of  the  strong  form  of  realizability  either 
removes  the  degrees  of  freedom  necessary  to  properly  calibrate  the  model  or  leads  to 
stiff  behavior. 

Despite  the  deficiencies  discussed  above,  the  SSG  model  seems  to  perform  moderately 
well  in  a  variety  of  homogeneous  turbulent  flows  as  documented  in  this  paper.  While  further 
improvements  would  be  most  welcome,  we  feel  that  it  is  unlikely  that  they  will  come  from  the 
standard  hierarchy  of  models  (61).  Fundamentally  new  approaches  will  be  needed.  Future 
research  will  be  directed  on  two  fronts.  The  SSG  model  will  be  implemented  in  a  full 
second-order  closure  for  the  computation  of  a  variety  of  complex  aerodynamic  flows  that  are 
of  technological  interest.  We  believe  that  when  the  SSG  model  is  used  within  the  framework 
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of  a  sound  second-order  closure,  it  may  be  possible  to  obtain  acceptable  engineering  answers 
for  a  range  of  turbulent  shear  flows  with  streamline  curvature.  In  parallel  with  this  effort, 
entirely  new  directions  in  modeling  the  pressure-strain  correlation  will  be  pursued.  These 
will  involve  the  introduction  of  a  tensor  length  scale  -  to  better  account  for  anisotropies  - 
and  the  possible  solution  of  a  transport  equation  for  Mx]kt  to  account  for  history  effects  in 
the  rapid  pressure-strain.  A  closer  examination  of  these  issues  will  be  the  subject  of  a  future 
paper. 
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APPENDIX  A 


Consider  the  tensor  function 

Pa  =  «y(b,S',^).  (41) 

Form  invariance  under  coordinate  transformations  (28)  requires  that  Fij  be  of  the  form  (see 
Smith  1971): 

F  =  Q!oI  4"  oqb  -b  o^b2  4"  0:3s  4"  O4S 

4-a5(bS'  +  S'b)  +  a6(b2S'  +  S'b2) 

4-a7(bS'2  +  S'2b)  4-  as(b2S'2  +  S'V) 

4-a9(bu^,  —  uT^b)  4-  aio^buT*  (A2) 

4-au(b  ui  —  a;  b  )  4-  ^12^^  bu>  —  u>  bu>  ) 

-|-Q'i3(S  U)  —  UJ  S  )  +  CL\^U  S  cj 

4-Q!i5(S  iJ1  —ZJ'S  )  4-  Q:j6(u7/S  u^,2  —  Zji2S  ZJ1)  4-  anZT12 
where  depend  on  all  possible  invariants,  i.e., 

aj  =  a,(/7,  III,tr  Z02,tr  S>2,tr  S>3,tr  bs'.fr  b2S,J 

tr  bS'2,tr  b2S'2,tr  buT'2,tr  b2^2,^  uT'bw''2, 

(**) 

tr  S'tJ^2,  tr  S'2u02,  tr  uFS  u02,  tr  bS  uT',  tr  b2S  uT1, 

tr  b^2S',  tr  bS  V),  i  =  0, 1, 2, ....  17. 

Taking  the  linear  part  of  Fa  in  S[j  and  w'j  yields  the  equation 

FW  =  p0l  +  pib  +  p2b2 +p3s' +  04{bS' 

4-S'b)4-/35(b2S'  +  S'b2)  {Al) 

4~/3g( buJ^  —  uT'b)  4"  /37(b2ui  —  u>  b2) 

where 

0o  =  000 (II,  III)  +  001  (II,  III)tr{ b  •  S')  +  0O2{II,  III)tr{ b2  •  S')  (A5) 

0i  =  0io (II,  HI)  +  Mil,  III)tr(b  •  S')  +  012(II,  Ill)tr{b2  •  S')  (A6) 
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(.47) 

(-48) 


02  =  020 (II,  III )  +  02,(11,  III)tr(b  •  S')  +  022(1I,  lll)lr(h2  ■  S') 

0i  =  0,(11,  II  I),  *  =  3, 4,-..,  7. 

Then,  by  taking  the  deviatoric  part  of  (since  Iljj  is  traceless)  and  multiplying  by  e  we 
obtain  II, y: 

n«  =  e(Jf  ’  -  1^%)-  (49) 

Equation  (29)  results  when  (A4)  -  (A8)  are  substituted  into  (A9)  and  the  identities  in 
equation  (27)  are  made  use  of. 


27 


APPENDIX  B 


The  pressure-strain  models  of  Shih  and  Lumley  (1985)  and  Fu,  Launder,  and  Tselepidakis 
(1987)  take  the  following  form: 

Shih-Lumley  Model 

'4 


n»j  =  —fiebij  4-^-4-  8a5j  K~5ij 


(ffiu  ■ 

ai + emfcjfir 


"4*7 'jk  _  4"  Emki^lrn^J  ^  TmnS  mn  &ij 

+  (3  +  Tas)  Tlk  {dTj  +  6mifcfim) 

^  4"  ^TmnS rnn^ij 


6  _  4  _  _ 

+  -TmnSmnkii  4*  —  [ 7,4 jlc  4"  TjkW ik] 


(51) 


4- 


5  K 


Tile  f  4"  r 


fdUi 

4 ~Tjk  I  4" 


Tkl  —  ^TnTjn^E  lm 


where 

tt5  =  _JL(1  +  OWF),  F  =  1  +  911  +  27111  (52) 

=  ///  =  ^babjkbki  (53) 

0  =  2  +  iF{80.Un[l  4-  62.4(-//  +  2.3///)]}  (54) 
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Fu,  Launder ,  and  Tselepidakis  Model 


-  -C*eVF  2 bij  +  47  (bikbkj  -  -bmnbmn6ij^J 


-0.6 


( ^ 
T\k  1  “I"  €mkj 


+Ti*  f‘HTL  +  emJtifiml  - 


Kdxk 

&v. 


+1'2Tm"  0^  "  Ik  {2t‘'t‘-s« 


— Ttt 


^*fc  I  "f"  n 


/ftFi 

Tjfe  (  ern.lSlr 


(55) 


+r{16//(rifcWi*  + 

~  24  bmi  bnj  ( TmkW nk  +  rnkWmk)} 

where 

Cl  =  -607/,  7  =  0.6,  r  =  0.7.  (56) 
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Equilibrium 

Values 

LRR  Model 

SSG  Model 

Experiments 

(^ll)oo 

0.193 

0.204 

0.201 

(622)00 

-0.096 

-0.148 

-0.147 

(612)00 

-0.185 

-0.156 

-0.150 

(SK/e)  oo 

5.65 

5.98 

6.08 

Table  1.  Comparison  of  the  predictions  of  the  Launder,  Reece  and 
Rodi  (LRR)  model  and  the  Speziale,  Sarkar  and  Gatski  (SSG)  model 
with  the  experiments  of  Tavoularis  and  Corrsin  (1981)  on  homogene¬ 
ous  turbulent  shear  flow. 
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Time  evolution  of  the  turbulent  kinetic  energy  in  rotating  shear  flows. 
Comparison  of  the  predictions  of  the  SSG  model  and  LRR  model  with 
the  large-eddy  simulation  of  Bardina,  Ferziger,  and  Reynolds  (1983)  for 
£o/SK0  =  0.296:  (a)  U/S  =  0,  (b)  fi /S  =  0.25,  and  (c)  tl/S  =  0.5. 


Time  evolution  of  the  turbulent  kinetic  energy  in  rotating  shear  flows. 
Comparison  of  the  predictions  of  the  SSG  model  and  LRR  model  with 
the  large-eddy  simulation  of  Bardina,  Ferziger,  and  Reynolds  (1983)  for 
eo/SKo  =  0.296:  (a)  O/S  =  0,  (b)  O/S  =  0.25,  and  (c)  0/5  =  0.5. 
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Time  evolution  of  the  turbulent  kinetic  energy  in  rotating  shear  flows. 
Comparison  of  the  predictions  of  the  SSG  model  and  LRR  model  with 
the  large-eddy  simulation  of  Bardina,  Ferziger,  and  Reynolds  (1983)  for 
e0/SKo  =  0.296:  (a)  0/5  =  0,  (b)  0/5  =  0.25,  and  (c)  0/5  =  0.5. 
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Comparison  of  the  models  with  large-eddy  simulations  (LES)  of  rotating 
shear  flow  for  fi/S  =  0,0.25  and  0.5:  (a)  LES  of  Bardina,  Ferziger,  and 
Reynolds  (1983),  (b)  FLT  model,  and  (c)  SSG  model. 
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Time  evolution  of  the  turbulent  kinetic  energy  in  plane  strain  for  e0/ SK 0  - 
2.0.  Comparison  of  the  predictions  of  the  LRR  model  and  SSG  model  with 
the  direct  simulations  of  Lee  and  Reynolds  (1985). 
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Time  evolution  of  the  anisotropy  tensor  in  plane  strain  for  eQ/SK0-  2.0. 
Comparison  of  the  predictions  of  the  LRR  model  and  SSG  model  with  the 
direct  simulations  of  Lee  and  Reynolds  (1985). 


1,RR  Model 


SSG  Model 


Simulations  O 


Figure  8.  Time  evolution  of  the  turbulent  kinetic  energy  in  plane  strain  for  eo/SKo  = 

1.0.  Comparison  of  the  predictions  of  the  LRR  model  and  SSG  model  with 
the  direct  simulations  of  Lee  and  Reynolds  (1985). 
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Time  evolution  of  the  anisotropy  tensor  in  plane  strain  for  co/SKq  =  1.0. 
Comparison  of  the  predictions  of  the  LRR  model  and  SSG  model  with  the 
direct  simulations  of  Lee  and  Reynolds  (1985). 
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Figure  10.  Time  evolution  of  the  turbulent  kinetic  energy  in  the  axisymmetric  con¬ 

traction  for  e0/SK0  =  0.179.  Comparison  of  the  predictions  of  the  LRR 
model  and  SSG  model  with  the  direct  simulations  of  Lee  and  Reynolds 
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Time  evolution  of  the  anisotropy  tensor  in  the  axisymmetric  contraction 
for  e0/SK0  =  0.179.  Comparison  of  the  predictions  of  the  LRR  model  and 
SSG  model  with  the  direct  simulations  of  Lee  and  Reynolds  (1985). 
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Figure  12.  Time  evolution  of  the  turbulent  kinetic  energy  in  the  axisymmetric  expan¬ 

sion  for  Eq/SKq  =  2.45.  Comparison  of  the  predictions  of  the  LRR  model 
and  SSG  model  with  the  direct  simulations  of  Lee  and  Reynolds  (1985). 
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Figure  13.  Time  evolution  of  the  anisotropy  tensor  in  the  axisymmetric  expansion  for 

Co/SKo  =  2.45.  Comparison  of  the  predictions  of  the  LRR  model  and  SSG 
model  with  the  direct  simulations  of  Lee  and  Reynolds  (1985). 
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